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I describe the magnetic-field and current-density distributions generated by two-dimensional (2D) 
pancake vortices in infinite, semi-infinite, and fmite-thickness stacks of Josephson-decoupled super- 
conducting layers. Arrays of such vortices have been used to model the magnetic structure in highly 
anisotropic layered cuprate high-temperature superconductors. I show how the electromagnetic 
forces between pancake vortices can be calculatated, and I briefly discuss the effects of interlayer 
Josephson coupling. 
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I. INTRODUCTION 

Since this paper is intended for publication in a special 
Festschrift issue honoring Mikc Tinkham, I have been in- 
vited to include some personal reflections in the introduc- 
tion. I believe I first heard his name when I was a gradu- 
ate student in the early 1960s at the University of Illinois- 
Urbana, working on extensions of the BCS theoryi to 
include anisotropy of the superconducting energy gap. 2,3 
A paper by Ginsberg, Richards, and Tinkham^ had re- 
ported results on the far-infrared absorption in supercon- 
ducting lead, which showed a precursor hump in the real 
part of the complex conductivity, ai(u))/aj^. I tried to 
explain this feature in terms of gap anisotropy but was 
unsucccssful. 

Throughout subsequent years, I have followed Mike 
Tinkham's career with considerable interest. I have ad- 
mired his research style, which consistently has resulted 
in new and interesting experimental results and theo- 
retical interpretations that advance the theory. I also 
admire anyone who can write carefully prepared books, 
and I have found his books on superconductivity (in 
both editions^) to be particularly uscful. I have asked 
students beginning research with me to work diligently 
through these books to learn the fundamentals of super- 
conductivity. 

One of the topics that Mike Tinkham finds interest- 
ing is vortex physics, and since this has been one of my 
main research interests, I would like to focus here on one 
aspect: two-dimensional (2D) pancake vortices. This is 
a favorite subtopic of mine, partly because I coined the 
name and partly because my 1991 paper on this subjectf 
has been so well received by the superconductivity com- 
munity (over 600 citations to date). Incidentally, al- 
though I wanted to put "2D pancake vortex" in the title 
of this paper, the editors of Physical Review B forbid 
this but did allow me to use these words in the abstract 
and the rest of the paper. I first reported on my work 
on 2D pancake vortices at a Gordon Research Confer- 
ence chaired by Mike Tinkham in June 1989, but (as has 
too often been the case with me) I was slow to publish, 
and some of the key results were published in 1990 by 
Artemenko and Kruglov 8 and by Buzdin and Feinberg, 9 



I later discovered that the basic solution had even been 
published in 1979 by Efetov^S but his work unfortunately 
had gone largely unnoticed. 

This paper is organized as follows. In Sec. II, I cal- 
culate the properties of 2D pancake vortices in an infi- 
nite stack of Josephson-decoupled superconducting lay- 
ers, first by considering all the layers as being very thin 
and then by considering the layers above and below the 
pancake layer as a continuumii In Sec. III, I use the 
continuum approach to calculate the properties of 2D 
pancake vortices in a semi-infinite stack of Josephson- 
decoupled superconducting layers. In Sec. IV, I again use 
the continuum approach to calculate the properties of 2D 
pancake vortices in a finite stack of Josephson-decoupled 
superconducting layers^ first for arbitrary thickness and 
then for a thickness much less than the in-plane penetra- 
tion depth, whcre the results bear some similarities to 
those of Pearliiiii^ for vortices in thin films. In Sec. 
V, I show how to calculate the electromagnetic forces be- 
tween pancake vortices, and in Sec. VI, I discuss some 
consequences of Josephson coupling. I conclude with a 
brief summary in Sec. VII. 



II. PANCAKE VORTEX IN AN INFINITE 
STACK OF SUPERCONDUCTING LAYERS 

The chief motivation for my work that led to the 
idea of the 2D pancake vortex was the question of how 
to describe the vortex structure of highly anisotropic 
layered cuprate high-temperature superconductors, with 
Bi-2212 (Bi2Sr2CaCu208-<5) being the best-known ex- 
ample. Applying the anisotropic Ginzburg-Landau 
eguationa 16 . 17 ! 18 ! 19 ! 20 . 21 ! 22 ! 23 ! 24 ! 25 ! 26 ! 27 ! 28 ! 29 to this mate- 
rial, it could easily be seen that the calculated value of 
the coherence length £ c (the length scale describing spa- 
tial variation of the order parameter in the c direction 
perpendicular to the layers) was less than the center-to- 
center distance s between the Cu02 bilayers. Since the 
Ginzburg-Landau theory assumes that all the character- 
istic lengths of superconductivity are large by comparison 
with atomic length scales, this fact indicated that some 
other theory was needed to describe details of the vortex 



2 



structure in the most anisotropic high-T c superconduc- 
tors. 

The natural way to incorporate the existence of dis- 
crete layers was to make use of the Lawrence-Doniach 
theory 30 which treats the intralayer behavior using 
Ginzburg-Landau theory but interlayer coupling via the 
Josephson effecti^i In this theory the coherence length 
£ c plays no role when its value is less than s, and the 
penetration depth A c describing the length scale of the 
spatial variation of supercurrents parallel to the c direc- 
tion can be related to the maximum Josephson supercur- 
rent J via 32 A c = (c4>o/8n 2 s Jo) 1 ^ 2 in Gaussian units. 
The parameter usually used to characterize the degree 
of anisotropy is 7 = A c /A a b, where A a & is the penetra- 
tion depth describing the length scale of the spatial vari- 
ation of supercurrents parallel to the layers (neglecting 
the anisotropy between the a and b directions, i. e., assum- 
ing for simplicity that A a « A/, m A a &). For Bi-2212, the 
value of 7 is so large that it is difficult to measure^ 7 was 
found in Ref. 0to be larger than 150, but a more recent 
quantitative determination^ 5 . has yielded 7 = 640 ± 25. 

For such highly anisotropic materials, it seemed sensi- 
ble to me to take the limit 7 —> 00 (A c = 00 or Jq = 0) 
as the starting point to describe vortex structure. The 
essential idea was that in a model of identical supercon- 
ducting layers separated by insulating layers, one could 
solve for the magnetic field and current density gener- 
ated by a 2D pancake vortex in one of the superconduct- 
ing layers when the other layers contained no vortices 
but served only to screen the magnetic field generated 
by the pancake vortex. With this solution as a building 
block, one could then find the magnetic field produced 
by a stack of such pancake vortices, even if misaligned, 
by the process of linear superposition. This was basi- 
cally the approach I had used in developing the theory 
that quantitatively explains the coupling forces between 
misaligned vortices in just two layerSf 3 ^! the primary 
and secondary superconducting layers of the dc trans- 
former studied experimentally first by Giaever; 3 ^ 3 ^ and 
later by Solomon^ Sherrill^i Deltour and Tinkhamj 4 ^ 
and Cladis et al.; 43 i 44 but in greatest detail by Ekin et 
ali 5 -^ 



A. Model of very thin discrete superconducting 
layers 

To calculate the magnetic-field and current-density dis- 
tributions generated by a pancake vortex in an infinitc 
stack of Josephson-decoupled superconducting layers, in 
Ref. I used the model in which the superconducting 
layers, all of thickness d, are centered on the plancs 
z — z n — ns (n = 0, ±1, ±2, ...), as sketched in Fig. 1. 
The London penetration depth within each layer is A s , 
such that the average penetration depth for currents par- 
allel to the layers is 3 ^ A|| = \ s (s/d) 1 / 2 , which corresponds 
to the penetration depth X a i, in the high-temperature su- 
perconductors. When the central layer (z = 0) contains 
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FIG. 1: Infinite stack of thin superconducting layers with a 
pancake vortex at the origin (bold arrow). 



a vortex at the origin but all other layers are vortex-free, 
the London fluxoid quantization condition 47 in layer n 
can be expressed as 

2>Kp[a<t > (p, z n ) + (2irA s /c)K < j,{p, z n )] = 4>oS n0 , (1) 

where in cylindrical coordinates a(p, z) = (fia^p, z) is the 
vector potential, K(p,z n ) = (j>K$(p,z n ) = $j^(p,z n )s is 
the sheet-current density in layer n averaged over the pe- 
riodicity length s, A s — 2A 2 |/s = 2A 2 /d is the 2D screen- 
ing length, and 4>q — hc/2e is the superconducting flux 
quantum. Equation QJ inevitably leads to a description 
of vortices in the London model- which is characterized 
by unphysical current-density and magnetic-field singu- 
larities on the vortex axis. The pioneering work on vor- 
tices by AbrikosoA* 4 ^ showed that such singularities are 
cut off at a distance of the order of the in-plane coher- 
ence length £ a f,. A simple model for the vortex core, em- 
ploying a variational core-radius parameter £ v ~ ^ a f,, has 
been used to describe straight vortices in isotropio^ and 
anisotropic 51 superconductors, as well as in films of arbi- 
trary thickness, whether isolatedS or in superconducting 
dc transformerss 3 - 7 » This model also could be used to cure 
the vortex-core singularities that are present in all the 
following results of this paper. 

If one takes the thickness d of each layer to be very 
small, as in Ref. 0, the vector potential can be expressed 
in the form 



a<t>(p, z) 



dqA(q)J 1 (qp)Z(q,z) 



(2) 



where J\(qp) is a Bessel function and Z (g, z) has scallops 
as a function of z that are necessary to describe the dis- 
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continuities of b p (p, z) arising from the induced sheet cur- 
rents K,p(p, z n ) for n^O. Note that b(p, z) = Vxa(p, z), 
such that 



and 



b p (p,z) 



b z (p,z) 



da^p, z) 
dz 

1 d\pa^{p,z)] 
p dp 



(3) 



(4) 



Inserting the exact expression for A(q) into Eq. yields 
a complicated integral that cannot be integrated analyti- 
cally. However, a close approximation to the exact result 
can be obtained by writing Z(q, z) = exp(—Q\z\), where 
Q = (q 2 + A^ 2 ) 1/2 and A(q) = (f) /2irA s Q; this approxi- 
mation, which is valid for s <C Aii, corresponds to retain- 
ing information on the scale of A|| but giving up detailed 
information on the finer scale of s. The resulting vector 
potential and magnetic field components are 



a<t>{p,z) 



b p (p,z) 



0oA| | 
27rA s p 



( e -|»IM|| _ e -r/A||) 5 



2ttA 



5_ e -r/A N 



&0 r± p -klA|| _ £„-r/A,| 



2nA s p l \z\ 



(5) 
(6) 
, (7) 



where r = (p 2 + z 2 ) 1 / 2 . Since in the high-temperature 
superconductors s/2Am = A||/A s ~ 10~ 2 , the vector po- 
tential term in Eq. JJJ of order Am/A s can be neglected 
in the central layer (n — 0), and we obtain to good ap- 
proximation 



K^(p,z ) 



ctpo 



47r 2 A s p 

However, for all the other layers (n ^ 0) we obtain 

C</> A|| 



K^p, Z n ) 



Air 2 A 2 p 



(e 



-i»/A, 



(8) 



(9) 



where z n = ns and r n = (p 2 + z 2 ) 1 / 2 . Note that the mag- 
nitude of the sheet-current density in the n = central 
layer is much larger, by a factor of order 10 2 , than the 
sheet-current density in one of the n ^ layers. It is for 
this reason that I gave the name pancake vortex to this 
field and current distribution. 

An interesting property of the above solutions is that 
the pancake- vortex-generated magnetic flux (fr z (p,z) = 
2irpa ( f > (p, z) up through a circle of radius p at height z is 
(using A s = 2A 2 |/s) 

$ z (p, z) = <Ms/2A||)( e -l*l/ A n - e- r / A N), (10) 

such that the magnetic flux up through a layer at height 
z is 



$*(oo,z) = </>o(s/2A||) e -l z l/ A n, 



(H) 



and the magnetic flux up through the central layer at 
z = is 



$*(oo,0) = 0o(s/2Ai 



(12) 



When s <C An, as in the high-temperature superconduc- 
tors, we see that <fr z (oo,0) 4>q. This is at first sur- 
prising until one realizes that fluxoids are quantized in 
superconductors but flux is not4i In the present prob- 
lem the fluxoid is the quantity on the left-hand side of Eq. 
JTJ, and since (2TrA s /c)K r f > (p, z ) is proportional to l/p 
and a^{p, zq) is very small, the fluxoid is due almost en- 
tirely to the current term. Note also that $^(oo, oo) = 0; 
this occurs because all the magnetic flux up through the 
central layer z = is directed radially outward by the 
screening currents in the layers with z > 0. 

On the other hand, an infinite stack of pancake 
vortices, whether straight or not, has quite different 
magnetic-flux properties. If there is one pancake vor- 
tex in every layer at z = z n = ns (n = 0, ±1, ±2, ...), 
then the magnetic flux up through the central layer (and 
by symmetry any other layer) is 



$ 2 (oo,O)=0 o (s/2A||) J2 



-M/A|| _ 



(13) 



where the last equality is obtained by evaluating the sum 
and making use of the property that s -C A||. Similarly, 
the radial magnetic field at p = oo and z = is now zero, 
since the positive contributions from all the pancake vor- 
tices below the central layer are canceled by the negative 
contributions from the pancake vortices above this layer. 



B. Continuum model 

The solutions given in Eqs. ©-© can be obtained 
more easily by regarding the n / layers as a contin- 
uum, characterized by the penetration depth A|| for cur- 
rents parallel to the layers^ Moreover, for a realistic 
treatment of stacks of just a few superconducting layers, 
a model accounting for finite layer thickness s is nccded. 
We therefore use the model sketched in Fig. 2 and write 
the London equationi£ in cylindrical coordinates in the 
form 

2Trp[a cj) (p, z) + (47rA 2 |/c)j0(p, z)} = (j) 5 n0 , (14) 

where the delta function on the right-hand side accounts 
for the presence of a vortex aligned along the z axis in 
the n = layer (|z| < s/2). Combining this equation 
with Ampere's law, — {c/ATr)(db p /dz — db z /dp) 7 and 
making use of Eqs. © and 10}, we obtain the partial 
differential cquation 



if- 



dz 2 dp 2 



1 da$ 
p dp 



1 1 



W = ~ ,2 S n0t (15) 
l'KAtp 
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FIG. 2: Continuum model of an infinite stack of supercon- 
ducting layers with a pancake vortex (bold arrow at the ori- 
gin) in the layer at z = zo — 0. 



which can be solved by writing a$(p, z) in the three re- 
gions z > s/2, —s/2 < z < s/2, and z < —s/2 in terms 
of Hankel componentsS as follows: 



a$(p,z) 



dqA a (q)J 1 (qp)e- Qz ,z>s/2, (16) 



o 



+ A -(q)e- 



o 2^A 2 Q 2 
+ Ao + {q)e Qz ] .h{qp), -s/2 < z < s/2, (17) 

dqA b (q)J 1 (qp)e- z ,z < -s/2, (18) 



where Q = (q 2 + l/A?.) 1 ' 2 . The four unknown functions 
A a (q), Ao-(q), Ao + (q), and Ab(q) can be obtained by ap- 
plying the boundary conditions of continuity of a<p(p, z) 
and b p (p,z) [Eq. ©] at the two interfaces z = ±s/2, 
carrying out the Hankel transforms using^i 



dppJMMtfp) = (l/q)5(q - q'), (19) 



and solving the four resulting linear equations. The re- 
sults are 



A _(q) = A 0+ (q) 



A a {q) = A b (q) 



4>o sinh(Qs/2) 
2ttA 2 Q 2 ' 



(20) 
(21) 



Note that s <C X\\, such that if we confine our attention 
to values of p ^> s, the integrals in Eqs. (|16|l - i|18|) are 



dominated by values of q <C l/s. We then may make 
the replacement sinh(Qs/2) — > Qs/2, which makes A a — 
Ab = 0o/27rA s <3, the same as A(q) in Ref. and Sec. II 
A. 

The magnetic flux up through a cirele of radius p in the 
plane with coordinate z is <& z {p, z) — 2irpa ( p(p, z). Evalu- 
ating the integrals for a${p, z) [Eqs. I|lt)|l - (|18|l ] in the limit 
as p — > oo, we can show without making the approxima- 
tion that s <C A|| that the pancake- vortex-generated mag- 
netic nux through a layer at height z, where \z\ > s/2, is 



$ z (oo,z) = o sinh(s/2A||)e~ |z|/A ii 



(22) 



and the total magnetic flux up through the central layer 
at z = is 

$ z (oo,0) = 0o(l-e- s / 2A n) «0o(s/2A,|). (23) 

If therc is one pancake vortex in every layer at z — 
z n = ns (n — 0, ±1, ±2, ...), even if they are misaligned, 
then by summing the contributions given in Eqs. 11' I'D 
and (|23H we find that the magnetic flux up through the 
central layer (and by symmetry any other layer) is exactly 
cf>Q. If all the vortices are aligned along the z axis, the 
magnetic-field and current-density distributions reduce 
to those of the London model, 48 for which 



bz(p) 



2np 

0o 



2np 

p 



1 - 



p 



K, 



P 



8vrA; 



P_ 

'An 



(24) 
(25) 
(26) 



and b p (p) = 0, where K n {x) is a modified Bessel funetion. 



III. PANCAKE VORTEX IN A SEMI-INFINITE 
STACK OF SUPERCONDUCTING LAYERS 

In Sec. II I reviewed the results found in Ref. for 
a pancake vortex in an infinite stack of supercondueting 
layers, where it is seen that the fields and currents decay 
exponentially on the scale of A|| above and below the 
layer containing the pancake vortex. For a sample of 
thickncss D ^> A| it is therefore elear that the fields and 
currents generated by pancake vortices that are many A|| 
from either surface are essentially the same as in Sec. II. 
However, the fields and currents are significantly altered 
when a pancake vortex is less than An from the surface 
of a sample of thickness D ^> A|| or when the sample 
thickness D is comparable with or smaller than An. In 
this seetion I use the continuum approximation deseribed 
in Sec. II B to obtain solutions deseribing the field and 
currents generated by a vortex in an arbitrary layer of a 
semi-infinite stack of supercondueting layers. In the next 
seetion (Sec. IV) I present solutions for a pancake vortex 
in a finite stack of arbitrary thickness D. 
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FIG. 3: Continuum model of a semi-infinite stack of super- 
conducting layers in the space z < with a pancake vortex 
(bold arrow) in the layer at z = zq. 



Consider a semi-infinite stack of superconducting lay- 
ers, with the top surface on the xy plane, such that all 
the layers are in the region z < 0, as sketched in Fig 
3. We number the superconducting layers such that the 
layer n = in the region zq- < z < zq + , centered at 
z = zq < 0, is the one containing the pancake vortex. 
The other layers are centered at z = z n — zo + ns, where 
positive (negative) n labels layers above (below) the pan- 
cake vortex. If there are N + layers above the pancake 
vortex, then the top layer is centered at z = zq + N + s — 
D/2 — s/2. By so numbering the layers, we still can use 
Eq. I|14|) as the London fluxoid quantization condition. 

As in Sec. II, we may write the vector potential in 
cylindrical coordinates as a{p,z) = (fra$(p,z). However, 
we now have different expressions for a<f,(p, z) in four re- 
gions: 



a<f,(p,z) = 

a<p(p,z) = 

a$(p,z) = 

a (p,z) = 



dqA > (q)J 1 (qp)e-i z ,z>0 1 (27) 

+A a+ {q)e^-^]J l {qp), 
zq+ < z < 0, 



+A Q+ ( q )e^ z -^}J 1 ( q p), 

z<3- < z < z 0+ , (29) 

dqA h (q)J 1 (qp)e Q{z - Zo \z< z -,(30) 



where Q = (q 2 + A, 2 ) 1/2 and z ± = z ± s/2. The six 
functions A>(q), A a -{q), A a+ (q), A -(q), A 0+ (q), and 
Ai,(q), obtained by applying the six boundary conditions 
of continuity of a$(ft, z) and b p (p, z) [calculated from Eq. 
©] at z — 0, zq + , and zq~, are 



A>(q) 
A a -(q) 
A a +(q) 
A -(q) 
Ao+(q) 

A b (q) 



0o sinh(Qs/2) 
ir\2Q 2 {l + q/Qy ' 

0o sinh(Qs/2) 
2ttA2Q 2 ' 

0o sinh(g£/2) , 1 -g/Q ^ 2 q z 



2ttA 2 Q 2 



4ttA 2 Q 2 



V 14 

Qs/2 

-Qs/2 
1 



q/Q' 



4ttA 2 Q 2 
-2sinh(Q.s/2)(i T ^) e 2 ^; 

0o sinh(Qs/2) 
2ttA 2 Q 2 



[1 



A-q/Q , 
A + q/Q' 



(31) 
(32) 

(33) 
(34) 

, (35) 
2Qz °].(36) 



Although the resulting integrals for a^(p, z) and thosc 
[via Eqs. © and J3J] for b p (p,z) and b z (p,z) can easily 
be calculated numerically, they are too complicated to 
evaluate analytically for arbitrary p and z. On the other 
hand, we can evaluate them approximately for large p. 
When p 3> Aji, the values of q that dominate the inte- 
grals in Eqs. (|2T|l - i|37)|l via the Bessel function Ji(qp) are 
those of order l/p <C l/Ay, such that we may replace all 
quantities under the integral except J\(qp) by their val- 
ues at q = 0. Similarly, because of the factor exp(— qz) 
in Eq. l(Po)l we may replace A>(q) by ^4>(0) to evaluate 
a^(p, z) when p is small but z ^> An. 

The magnetic flux up through a circle of radius p in 
the plane with coordinate z is & z (p,z) — 2 r Kpa r f,(p,z). 
Evaluating the integrals as indicated above for a^p^z) 
in the limit as p — > oo, we obtain for the total magnetic 
flux up through the plane with coordinate z: 



$ 2 (oo,z) = 20 o sinh(s/2A||)e 2o/A n ;Z > o, (37) 
$ z (oo,z) = 20 o sinh(s/2A||)cosh(z/A||)e Zo/A n ; 

zq+ <z < 0, (38) 
$,(oo,z) - o {l-cosh[(z-z o )/A||] e - s / 2A M 
+ sinh(s/2A||)e( z+ ^/ A M}, 
(28) z o- < z < z Q+ , (39) 

$> z {oo,z) = 20 o sinh(s/2A||)cosh(z o /A||)e z/A ii, 

z<z _. (40) 



When the pancake vortex is in the top layer (i. e., when 
zq = —s/2), the magnetic flux $ z (oo,0) up through 
the top surface is approximately 0o(s/A||), since s/Ah ~ 
10~ 2 <C 1. When the pancake vortex is in a layer much 
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farther than A|| from the top surface, the amount of mag- 
netic flux up through the top surface $ z (oo, 0) [Eq. (|3*7|) ] 
becomes exponentially small (recall that zq < 0). The 
precise magnetic ficld distribution generated in the space 
above the superconductor within A| | of the origin can be 
calculated numerically for a given pancake-vortex posi- 
tion z from Eqs. 10, |@J|, and (|T7Jl . However, at dis- 
tances r = \J 'p 2 + z 2 somewhat larger than Ai from the 
origin, we have to good approximation for z > 



b P (p,z) 
bz(p,z) 



$*(oo>0) (1 _ z 
2np r 

$ z (oc,0) p 
2n r 3 ' 

$z(°°,0) z 
2n r 3 



(41) 
(42) 
(43) 



In other words, the magnetic field generated by the 
pancake vortex appears as if generated by a magnetic 
monopole, with the fhrx <I> 2 (oo, 0) [Eq. (|37[) ] spreading out 
evenly into the hemisphere above the surface. It is im- 
portant to note that only pancake vortices within about 
A|| (or A a b in the high-temperature superconductors) 
are visible using Bitter decoration, scanning Hall-probe 
microscopy, scanning SQUID microscopy, or magneto- 
optical techniques; pancake vortices deeper than this 
make an exponentially small contribution to the mag- 
netic field above the surface. 

From Eq. 1|39|1 we see that the magnetic fiux up through 
the plane z = zq in the layer containing the pancake 
vortex is 

$ 2 (oc,z ) = 0o[l - e - s/2A n + sinli(s/2A||)e 2 " :o/A ii] 

w ^( S /2A||)(l + e 2z °/ A ii). (44) 

When the pancake vortex is in the top layer (i. e., if 
z = —s/2), the magnetic flux up through this layer is 
approximately 4>o(s/\\\), and when the pancake vortex is 
deep inside the superconductor (i. e., if — zq S> A||), the 
magnetic flux up through the pancake layer is approx- 
imately </>o(s/2A||), as found in Sec. II for the infmite 
superconductor [Eqs. I|12|l and 1)23(1] . 

If there is a pancake vortex in every layer, even if they 
are misaligned, the total magnetic flux up through any 
plane with coordinate z is exactly equal to 0o- This can 
be shown by replacing zq by z n = zq + us, noting that the 
top layer is centered at —s/ 2, and summing over all n, 
using Eq. I|37|) if z > 0. On the other hand, if z < 0, one 
must use Eq. I|4(J|) for the top layers for which z n — s/2 > 
z, Eq. 1|39[) for the layer containing z for which z n — s/2 < 
z < z n + s/2, and Eq. I|38|l for the remaining layers for 
which z n + s/2 < z. If all the pancake vortices are aligned 
along the z axis, the magnetic-field and current-density 
distributions reduce to those calculated by PearliiiS 3 ^ 
for a vortex in a semi-infinite superconductor. 

Scanning Hall-probe experiments visualizing vor- 
tices in underdoped, highly anisotropic YF^CuaOg+z 
(YBCO) single crystals, where x — 0.35 — 0.375, re- 
cently have been carried out by Guikema. 5 £ In the most 



underdoped crystals, the observations revealed what at 
first appeared to be "partial vortices" carrying magnetic 
flux less than 0o- Guikema concluded, however, that 
such images are caused by a full vortex that is par- 
tially displaced horizontally, i. e., a "split pancake- vortex 
stack." The magnetic flux generated above the surface 
by the two parts of the vortex stack can be calculated 
as follows. Suppose the bottom portion, consisting of 
pancake vortices below the plane z = —d, is aligned 
along the z axis, and the top portion, consisting of pan- 
cake vortices above the plane z — —d, is aligned par- 
allel to the z axis but at (x,y) = (xq,0). Using Eq. 
(|37|l to sum the contributions from the pancake vortices 
in the two portions, one finds that the magnetic flux 
&bot = 0o exp ( — d/ A 1 1 ) generated by the bottom portion 
emerges from the vicinity of the origin (x, y, z) = (0, 0, 0), 
and the magnetic flux $t op = 4>o[l — exp(— <i/A||)] gener- 
ated by the top portion emerges from the vicinity of the 
point (x,y,z) — (xo,0, 0). The two fiux contributions 
should be resolvable when the displacement xq exceeds 
the Hall-probe size and the probe's field sensitivity allows 
detection of both contributions. 



IV. PANCAKE VORTEX IN A FINITE STACK 
OF SUPERCONDUCTING LAYERS 



Since all laboratory samples of the high-temperature 
superconductors are of finite thickness, it is important to 
examine how the properties of pancake vortices discussed 
in Sec. II are modified when we take the finite thickness 
into account, including the possibility that the thickness 
D may be less than the penetration Am. Let us begin by 
considering a pancake vortex centered on the z axis in a 
stack of superconducting layers, each of thickness s, in 
the region —D/2 < z < D/2, as sketched in Fig. 4. As in 
the previous section, we number the layers such that the 
layer n = at z = Zq, where \zo\ < (D/2 — s/2), is the 
one containing the pancake vortex. The other layers are 
centered at z = z n = zq + ns, where positive (negative) n 
labels layers above (below) the pancake vortex. If there 
are N + layers above the pancake vortex, then the top 
layer is centered at z = z + N + s = D /2 — s/2, and if 
there are JV_ layers below the pancake vortex, then the 
bottom layer is centered at z = zq — N-S = —D/2 + s/2. 
As in Sec. II B, I treat all the layers using the continuum 
approximation and use Eq. (|14J) as the London fluxoid 
quantization condition. In Sec. IV A, I show how to 
caleulate the fields generated by a pancake vortex in a 
finite stack of Josephson-decoupled superconducting lay- 
ers, each of thickness s, with an arbitrary total stack 
thickness D relative to A||. In Sec. IV B, I consider the 
simplifications that arise when D <^ An, which corre- 
sponds to the case of high-temperature superconducting 
samples consisting of roughly ten or fewer unit cells along 
the c direetion. 
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FIG. 4: Continuum model of a stack of superconducting layers 
in the space \z\ < D/2 with a pancake vortex (bold arrow) in 
the layer at z — zq. 



A. Finite stack of arbitrary thickness 

As in Sees. II B and III, we write the veetor potential 
in cylindrical coordinates as o(p, z) = <pa^,(p,z). How- 
ever, we now have different expressions for a$(p,z) in 
five regions: 



a<f>(p,z) = 

a<t>(p,z) = 

a<t>{p,z) = 

a<p(p,z) = 

a (p,z) = 



dqA > {q)J 1 {qp)e-^- D l' 1 \ 
z > D/2, 

dq[A^( q )e-^- D W 



(45) 



+A a+ { q )e^ z - D ^]J 1 (qp) 1 
z + < z < D/2, 



(46) 



+ A -(q)e-^ z - z ^ 



+A 0+ (q)e^-^]J 1 (qp), 

zo- < z < z 0+ , (47) 

dq[A h _(q)e- < K* +D M 

+A b+ (q)e^ z + D ^}J 1 (qp), 
-D/2<z< zo-, (48) 

/ dqA < {q)J 1 {qp)e^ z + D l 2 \ 

z < -D/2, (49) 



where Q — (q 2 +A, 2 ) 1 ^ 2 , the subseript a (b) denotes the 
layered region above (below) the pancake vortex, and 



z ± = z ± s/2. The eight funetions A > (q), A a _(q), 
A a+ (q), A _(q), A 0+ (q), A b _{q), A b+ (q), and A<(g), 
obtained by applying the eight boundary conditions of 
continuity of a ( p(p, z) and b p {p,z) [caleulated from Eq. 
©] at z = D /2, z -, z 0+ , and —D /2, are 



A>{q) 
A a -{q) 
A a +{q) 
A -(q) 



n\\Q 2 



00 



2n\2Q 2 



2n\ 2 Q 2 



smh(Qs/2)G(q,z ), (50) 



sinh(Q S /2)(l + q/Q)G( q ,z ), (51) 



47rA 2 Q 21 



smh(Q S /2)(l-q/Q)G( q ,z ), (52) 



2 sinh(Qs/2)(l + q/Q)G(q, zo)e QD/2 e~ Q2 °] 



00 



47rA 2 Q 2 ' 



-Qs/2 



2smh(Qs/2){l - q /Q)G(q,-z )e-^ D / 2 e~ Qzo }, 

(53) 



47rA 2 Q 2 



r e -Qs/2 



2 sinh(Qs/2) (1 - q/Q)G(q, z )e 



-QD/2 e Q Zo i 



47rA 2 Q 2 



2 f e Qs/2 



A+(<z) 

A<(q) 
where 
G(q,z) 



2 sinh(Qs/2)(l + q/Q)G(q, -z )e QD / 2 e Qz °] 

(54) 

swh{Qs/2)(l-q/Q)G(q,-z ), (55) 



27rA 2 Q 2 



27rA 2 Q 2 



smh(Qs/2)(l + q/Q)G{ q ,-z ), (56) 



ttA 2 Q 2 



sinh(Qs/2)G(q, -z ), 



(57) 



(l + q/Q) + (l-g/Q)e-Q D e- 2 Q z 
(l + q/Q) 2 -(l-q/Q) 2 e- 2 Q D 



e -QD/2 e Qz_ 



(58) 

Although the resulting integrals for a^(p, z) and those 
[via Eqs. © and |@J] for b p (p,z) and b z (p,z) can easily 
be caleulated numerically, they are too complicated to 
evaluate analytically for arbitrary p and z. On the other 
hand, we can evaluate them approximately for large p. 
When p 3> Am, the values of q that dominate the inte- 
grals in Eqs. (|45l) - 149fl via the Bessel funetion J\{qp) are 
those of order l/p <C 1/Au, such that we may replace all 
quantities under the integral except J\{qp) by their val- 
ues at q — 0. Similarly, because of the faetors exp(— qz) 
and exp(gz) in Eqs. I|45|l and (|49|l we may replace A > (q) 
by ^4>(0) and A < (q) by A<(0) to evaluate a<f,(p, z) when 
p is small but \z\ — D/2 ^>> A|| . 
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The magnetic flux iip through a circle of radius p in in the limit as p — > oo, we obtain for the total magnetic 
the plane with coordinate z is & z (p,z) = 2irpa ( / ) (p, z). flux up through the plane with coordinate z>— 
Evaluating the integrals as indicated above for a<f,(p,z) 



$ z (oo,z) = 24> sinh (-i-) cosh ( D ' 2 + *° ) / sinh (^-),z>D/2, (59) 



' 2A|| Aii Air 



$ z (oc,z) = 20 o sinh (-f-) cosh ( D, \ + Z ° ) cosh Z ) / sinh A, z 0+ < z < D /2, (60) 



'2A|| A|| A|| 

$ z (oo, z) = O {1 - [e^- s )/ 2A n cosh (£I1ZH) _ e -(D-s)/2X n CQsh ( g/2±gD ) ]e »An /2 sinh (£_ 

A,, A,, A,, 

_ [e (D-s)/2X n CQsh ( J>/2 + Z ) _ g-CD-J/aA,, CQsh ( J>/2 -Z0 )]e _ zA|/2s . nh ( ^ }) 

A|| A|| Aj| 

zo- < ^ < 20+, (61) 

$ z (oo,z) = 20 o sinh(-^-)cosh( i: '^ 2 " 2:o )cosh( £> ^ + Z )/sinh(^-),-£'/2 < z < z -, (62) 

II II II II 

$ z (oo,z) = 20 o sinh (-£-) cosh ( D, \ ~ Z ° )/ sinh A , z < -D /2. (63) 

1 1 II II 



The magnetic flux $ z (oo,Z3/2) up through the top 
surface is given by Eq. Whcn Z) ^> A|| and a 

pancake vortex is in the top layer (i. e., when zq = 
D/2 - s/2), we obtain $ z (oo,D/2) w o (s/A||), which 
is a tiny fraction of 0o, since s/A| ~ 1CP 2 -Cl. As a 
function of the distance D/2 — zq of the pancake vor- 
tex from the top surface, we see that $ z (oo,D/2) ~ 
0o(s/A||)exp[-(£)/2 - z )/X\\]. When D < A||, we find 
that <& z (oo, D/2) w <fi (s/D) = (f>o/N, independent of the 
position zo of the pancake vortex within the stack, where 
N = D/ s is the number of layers in the sample. When 
N = D/s = 1, $ z (oo,Z?/2) = 0o, because our results 
then reduce to those of Pearljiiiii^ who calculated the 
field and current distribution generated by a vortex in a 
film of thickness much less than the London penetration 
depth. The precise magnetic field distribution generated 
in the space above the superconductor can be calculated 
numerically for a given zo from Eqs. |JS}, (@J, and (|45ll . 
However, at distances r + = y/ p 2 + (z — D/2) 2 from the 
point on the surface directly above the pancake vortex 
that are larger than A|| when D > 2A|| or larger than 
the two-dimensional screening length A d = 2A 2 / D when 
D < 2A||, we have to good approximation for z > D /2 



a<j,(p,z) 
b p (p,z) 
bz(p,z) 



*.(oo, P/2) (z -D/2) 

2tt P 1 r+ J ' v ' 

$ z (oo,£>/2) p 

2tt r\> 
<$> z (oo,D/2) jz~D/2) 

2tt r? 



(65) 
(66) 



cake vortex appears as if generated by a positive magnetic 
monopole, with the flux $ z (oo, D/2) [Eq. (|59|) ] spreading 
out into the hemisphere above the surface. 

Similar statements can be made about the magnetic 
flux <E> z (oo, — D/2) up through the bottom surface [Eq. 
|JH3}]. At large distances r_ = y/ p 2 + (z + D/2) 2 from 
the point on the surface directly below the pancake vor- 
tex, the magnetic field appears as if generated by a neg- 
ative magnetic monopole. 

From Eq. (|61[) we see that the magnetic flux up through 
the plane z = zo in the layer containing the pancake 
vortex is 



$ z (oo,z ) = o {l-[sinh( 



D -s/2 , 
A,, ' 



In other words, the magnetic field generated by the pan- 



- sinh (^T-) cosh (t^)]/ sinh (t)^ 

2A|| A|| A|| 

(67) 

When D ^> Am, the dependence of this magnetic flux 
upon the distance (D /2 — \zq\) from the top or bot- 
tom surface is given by $ z (oo,zo) ~ <j>o(s/2X\\){l + 
exp [— 2{D/2 — l^oD/An]}. When the pancake vortex is 
in the top or bottom layer (i. e., if \zq\ = D /2 — s/2), 
the magnetic flux up through this layer is approximately 
0o(s/A||), and when the pancake vortex is deep inside 
the superconductor (i. e., if D/2 — \zq\ ^S> X\\), the mag- 
netic flux up through the pancake layer is approximately 
0o(s/2A||), as found in Sec. II for the infinite supercon- 
ductor [Eqs. JT2J and (ESJ)]. When D < A||, we see that 
<j> z (oo,z ) « 4>o{s/D) = (f>o/N, independent of the po- 
sition zo of the pancake vortex within the stack, where 
N = D / s is the number of layers in the sample. 
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If there is a pancake vortex in every layer, even if they 
are misaligned, the total magnetic flux up through any 
plane with coordinate z is exactly equal to 0o- This can 
be shown by replacing zq by z n = zq + ns and summing 
over all n, using Eq. (ISU) if z > D/2 or Eq. (JB^J if z < 
—D/2. On the other hand, if \z\ < D/2, one must use 
Eq. (|62f> for the top layers for which z n — s/2 > z, Eq. 
(I61|l for the layer containing z for which z n — s/2 < z < 
z n + s/2, and Eq. (UU for the remaining layers for which 
z n + s/2 < z. If all the vortices are aligned along the z 
axis, the magnetic-field and current-density distributions 
reduce to those given in Ref.|E3 when = 0. 

It is possible that scanning Hall-probe or magneto- 
optical experiments may be able to detect partial vor- 
tices or split pancake-vortex stacks^ carrying magnetic 
fiux less than 0o in samples of highly anisotropic lay- 
ered superconductors of thickness D < Aii. The mag- 
netic flux generated above the surface z = D /2 by the 
two parts of the vortex stack can be calculated as fol- 
lows. Suppose the bottom portion, consisting of pan- 
cake vortices below the plane z = D /2 — d, is aligned 
along the z axis, and the top portion, consisting of pan- 
cake vortices above the plane z = D /2 — d, is aligned 
parallel to the z axis but at (x,y) = (xq,0). Using 
Eq. Ij59(l to sum the contributions from the pancake 
vortices in the two portions, one finds that the mag- 
netic ffux $ bot — o sinh[(Z? — d)/A||]/sinh(D/A||) gen- 
erated by the bottom portion emerges from the vicin- 
ity of the point (x,y,z) — (0,0, D/2), and the magnetic 
ffux $ top = O {1 - sinh[(D - rf)/A||]/sinh(Z»/A||)} gen- 
erated by the top portion emerges from the vicinity of 
the point (x,y,z) — (xq,0,D/2). The two flux contri- 
butions should be resolvable when the displacement xq 
exceeds the Hall-probe size and the probe's field sensi- 
tivity allows detection of both contributions. Note that 



*6ot = *to P = 0o/2 when d = D/2 < A 



II- 



B. Finite stack of thickness D <C Ai 

Considerable simplifications occur when the thickness 
D = Ns of the stack is much less than the in-plane pen- 
etration depth A||ii It is well known from the work of 
Refs. and IhH that when D < An the characteristic 



screening length in isolated films is not A|| but rather 
the 2D screening length Ar, = 2X^/D. This is also 
true for the case of Josephson-decoupled stacks of to- 
tal thickness D considered here. We may derive equa- 
tions for a^,(p,z), b p (p,z), and b z (p,z) valid for D <C A| 
and p » A|| by starting with Eqs. ijlB^ - lf^ . applying 
Eqs. J3J and (@J, and making the replacement e ±( ^ z = 
cosh(Qz) ± sinh(Qz). Since we are most interested in 
values of p of the order of or larger, because of the 
presence of J\(qp) the dominant values of q in the re- 
sulting integrals are of the order of q ~ 1/A^ <C 1/An, 
such that Q can be replaced by 1/An, and small quanti- 
ties of the order of D/A\\ and qX^ are of the same order 
of magnitude. Expanding in powers of the small quan- 
tities (D/A|| and <?A||), we find that both a<p{p, z) and 
b z (p,z) are to lowest order independent of z, with small 
correction terms of the order of D/X\\, such that to good 
approximation we may write these quantities as 



a<t>{p,z) 
a<t>(p, z) 
a<t>(p,z) 
bz{p,z) 
b z (p,z) 
bz{p,z) 



dq 



J\{<IP) -q(z-D/2) 



2ttN J 1 + qA D 



z > D/2, 



(68) 



< i dq JlM ,-D/2<z<D/2, 



2nN J q 1 + qA D 



" i dq JiML e ^+D/z) jZ <- D/2 , 



(69) 



2nN J q 1 + qA D 



(70) 



dq 



qJo(qp) 



-q(z-D/2) 



2ttN J l + qA 



>D/2, 



" i dq ^kA : _ D/2 < z < D/2 , 



2ttN J 1 + qA L 



2-kN J * 1 + qA D 



(71) 
') 

(72) 
12, 
(73) 



On the other hand, the radial component of the magnetic 
field varies strongly with z: 



b P {p,z) 


2irN J 


b p {p,z) 


00 r° 

2ttN J q 


b P {p, z) 


= b p (p,D/2) + 


b P {p,z) 


= b p (p,D/2) + 



gJl(gp) e -«.-p m >D/2 

gJl(g f } ,z = D/2, 
l + qA D ' ' ' 

(D/2 -z) a^p) <z<D/2 

-(D/2) X^' Z0+ - Z - D/2 > 

(D/2 - z) a^p) _ (z 0+ - z) 0Q 

(D/2) A D (D/2) 2kA dP 



zq- < z < Z Q+ , 



(74) 
(75) 
(76) 
(77) 
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b P {p,z) 
b P {p,z) 

b P {p,z) 



b p (p,D/2) 



(D/2 - z) a (p) _ s fa 
(D/2) A D (D/2) 2nA D p 



-D/2 < z < z _, 



b p (p,-D/2),z=-D/2, 

r^siM -p/2, 

2irNj Q q l + qA D ~ 1 



(78) 

(79) 
(80) 



where we use a^{p) to denote the vector potential in 
the region \z\ < D/2, since a^(p, z) is very nearly inde- 
pendent of z. The sheet current K n (p) = K$(p,z n ) — 
sjtf,(p, z n ) in layer n can be obtained from cithcr 
j#(p, z) — (c/ 'A-n)db p (p, z)/dz or the fhixoid quantization 
condition [Eq. iflijl]: 



K n {p) 



2ttA s l 2ttp 



[w^Z S nO ~ a <t>(p)] ■ 



(81) 



where A s = 2Am/s = NAd- The net sheet current 
through the thickness D is the sum of the K n : 



K D (p)= K ^P) 



-N— 



2txA d y 2-nNp 



i^W-~ a M- (82) 



The integrals appearing in Eqs. (|68[) - l|82[) . which are 
evaluated in Appendix A, have simple forms in the lim- 
its D <C p <C Ad and p ^> Ae>. The corresponding ex- 
pressions for the physical quantities we have calculated 
in this section are given in Table[I] The magnetic-field 
and current-density distributions reduce to the thin-film 
results of Pearli^iii when N — 1 and D — s or when each 
of the N layers contains a pancake vortex on the z axis. 



TABLE I: Results for one pancake vortex in a stack of N 
superconducting layers of total thickness D = N s -C An in 
the limits D < p -C A_d = 2A^/£> and p » A D . Sinco 
D is very small, r = (p 2 + z 2 ) 1 ^ 2 may be regarded as the 
distance from the pancake vortex, and \z\ may be regarded as 
the distance from the top or bottom surface. 



Physical quantity 

a<t>(p,z) 

a<t>(P,°) 
$ z (p, z) = 2npa ci> (p, z) 
$ z (p,0) = 2np a4 ,(p,0) 
b p (p,z), z = ±\z\ 
b p (p, ±D/2) 
b z (p,z) 
6«(p,0) 
Ko(p) = K<t,(p,zo) 
K n (p) = K$(p,z n ), n 
Kd(p) 



p < A D 

<i>o( r -\z\) 
2irNA D p 

4>a 
2ttNK d 
<t>a( r -\ z \) 
NA D 

<PoP 
NA D 
j_ <t> (l — \x\) 



" 2nNA D pr 



± 



4>a 



2irNA D p 

2ttNA d t 
<t>a 

2ttNA d p 

4tv 2 NA d p 

ct t>0 

4tv 2 N 2 A 2 d 

c<t>0 
4iv 2 NA D p 



P » A D 

2-KNpr 

4>0 
2-KNp 
4>o(r-\z\) 
Nr 

'i f) 

N 

-\-_iop 

2wNrX 

+ <t>a . 
2-irNp 2 

<t>QZ 



2wNrZ 
4>0>~D 
2-rrNp 3 
c4>q(N-1) 
4tv 2 N 2 A d p 

— c tP 
4tv 2 N' 2 A dp 

c<f>0 
4ir 2 Np 2 



semi-infinite stack of superconducting layers is given to 
excellent approximation by 



V. FORCES 

The force on a second pancake vortex at the position 
(p, z n ) due to a pancake vortex centered on the z axis at 
(0, Zq) can be calculated from the Lorentz force^ Since 
pancake vortices cannot move out of their planes, the 
force is directed parallel to the planes in the radial (p) 
direction: 



f p(p) = K^{p,z n )(j) /c, 



(83) 



where 



K$(p, z n ) 



C r (p0 x 

2ttA s 2np 



a<f,(p,z n )] (84) 



is the sheet-current density and a<jj(p, z n ) is the vector 
potential at (p, z n ) generated by the pancake vortex at 
(0, Zq), and A s = 2Aj 2 | /s = NA D . 

If both pancake vortices are in the same plane, the 
interaction force is always repulsive and in an infinite or 



Fp(p) 



An 2 A s p 



(85) 



for all p. The reason for this is that the vector potential 
in Eq. (jg) obeys a (p, z ) < (s/A||)(0 Q /27rp) < 4>o/2ttp, 
as shown in Sees. II and III. However, for a finite stack of 
thickness D <C Am consisting of N layers, Eq. I|85|l holds 
only for small p (p< Ac), where the vector potential in 
Eq. H84JI is much smaller than cf>o/2irp. As diseussed in 
Sec. IV A, the magnetic flux at infinite radius $2(00, zq) 
up through the pancake-vortex layer is approximately 
<po/N, which means that a^p, zq) w 4>q/2ttNp for large 
p, and 



F P {p) 



(N - 1) 0; 



N 4:ir 2 A s p 



,P> Ai 



(86) 



In the special case when N = 2, the repulsive force given 
in Eq. is half that in Eq. 

If the two pancake vortices are in different planes, the 
4>o/2irp term in Eq. (|84|) is absent, and the interaction 
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force is given by 



F P {p) 



4>oa<t,{p, z n ) 
2?rA, 



(87) 



Because a<p(p, z n ) is always positive, the interaction force 
is always negative, i. e., in a direction so as to cause the 
two pancake vortices to become aligned along the same 
vertical axis. For the general case, it is not a simplc 
matter to calculate the spatial dependence of the attrac- 
tive force between pancake vortices in different layers, 
as can be seen from the expressions for a^(p, z) given in 
previous sections. However, we can say that for an infi- 
nite or semi-infinitc stack of superconducting layers, the 
magnitude of this attractive force is orders of magnitude 
smaller than the repulsive force between pancake vortices 
in the same layer. The attractive force between vortices 
in different layers in an infinite stack (or deep inside a 
semi-infinite stack) has a range A|| in the z direction. 
Equation JSJ) shows that the attractive force in the in- 
finite stack vanishes exponentially when the interplanar 
separation of the pancakes along the z direction exceeds 
An. For a finite stack of thickness D — N s <C An, we 



find that the attractive force between pancake vortices 
in different layers is 



F P {p) 



4tt 2 A2 



An 2 N 2 A 2 D 



D -C p < A D , 



which agrees with the force in the infinite stack caleulated 
from Eq. J^J when \z n \ <C p <C Aii, and 



F P (p) 



4tt 2 N 2 A dP 



p> A D . 



(89) 



For the special case of two layers (N = 2) and a sep- 
aration p » Ao, the magnitude of the attractive force 
exerted by a pancake vortex in one layer upon a pancake 
vortex in the other layer [Eq. (|89|l ] is equal to the magni- 
tude of the repulsive force between two pancake vortices 
in the same layer [Eq. with A s = NAd]- 

The energy per unit length of a uniformly tilted infi- 
nite stack of pancake vortices in an infinite stack of su- 
perconducting layers was caleulated in Ref.0 The corre- 
sponding line tension T(0) was caleulated in Ref. as a 
funetion of the angle 9 relative to the z axis and found to 
be positive only for 9 < 51.8°, indicating an instability 
beginning at 51.8°. Further calculationsS showed that, 
because pancake vortices energetically prefer to line up 
parallel to the z axis, the energy for an infinite stack of 
pancake vortices with a large average tilt angle is reduced 
when the stack breaks up into shorter staeks parallel to 
the z axis with kinks between them. Pe et al^ caleulated 
the equilibrium positions of a stack of pancake vortices in 
a finite stack of Josephson-decoupled layers when equal 
and opposite transport currents are applied to the top 
and bottom layers. They found that the pancake vortices 
in the top and bottom layers have large displacements to 
the left and right, while the other vortices all remain 



elose to the z axis. Related model caleulations were car- 
ried out in Ref. ED for moving two-dimensional pancake 
vortex lattices in a finite stack of magnetically coupled 
superconducting thin films with transport current only 
in the top layer. For small currents, the entire electro- 
magnetically coupled vortex array moves uniformly in 
the direction of the Lorentz force but with a large dis- 
placement of the pancake vortices in the top layer rela- 
tive to the others, which remain in nearly straight lines 
perpendicular to the layers. Above a critical decoupling 
current, the 2D vortex array in the top layer periodically 
slips relative to the arrays in the other layers, and the 
dc current- voltage charaeteristies for the top and bottom 
layers resemble those reported in Refs. and EH for the 
dc transformer. 



VI. JOSEPHSON COUPLING 

The equations underlying the solutions presented in 
Sees. II, III, and IV assume no interlayer Josephson cou- 
pling. Implicit in these solutions is the assumption that 
the component of the magnetic field parallel to the lay- 
ers spreads out uniformly in the radial direction. This 
is consistent with the idea that if a magnetic field is ap- 
plied parallel to a stack of Josephson-decoupled layers, 
the field will penetrate uniformly between the layers. 

When the layers are Josephson-coupled, however, par- 
allel magnetic fields penetrate the structure in the form 
of quantized Josephson vorticesp^S As diseussed in Ref. 
Ifi2[ Josephson vortices in the high-temperature supercon- 
duetors have highly elliptical current and field patterns. 
Since the decay length for currents perpendicular to the 
layers is A c and that for currents parallel to the layers 
is A a b, the ratio of the width of the pattern parallel to 
the layers to the height perpendicular to the layers is 
7 = A c /A a (, at large distances from the non-linear Joseph- 
son core. For a high-K Abrikosov vorteji^ in an isotropic 
superconduetor, the decay length at large distances is 
the penetration depth A, and the currents in the non- 
linear core vary on the much smaller length scale of the 
coherence length £. The behavior in a Josephson vortex 
is analogous. The small length scale for spatial varia- 
tion of the Josephson currents in the vortex core (whose 
axis is centered in the insulating layer between two adja- 
cent superconducting layers) is the Josephson lengthSiS^ 
A,/ = js, while the corresponding length scale for the 
return of these currents parallel to the layers is s, such 
that the ratio of the width to the height of the Josephson 
core is 7 = Aj/s = A c /A ab . 

In the presence of interlayer Josephson coupling, the 
magnetic-field and current-density distributions gener- 
ated by a pancake vortex are unaltered at short distances 
but are strongly affected at distances of the order of A,/ 
and A c . To give an example, imagine an infinite stack of 
semi-infinite Josephson-coupled superconducting layers, 
all parallel to the xy plane, filling the half-space x > 0, 
such that the surface coincides with the plane x = 0. 
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Imagine creating a pancake vortex at the origin in the 
superconducting layer n — and moving it in to a dis- 
tance xq. The magnetic-field and current-density distri- 
butions, including the effects of a dipole-like stray field 
that leaks out into the space x < 0, have been calcu- 
lated as a function of xq in Ref. |3 under the assump- 
tion of very weak Josephson coupling. In the presence of 
Josephson coupling, however, the component of the mag- 
netic field parallel to the layers cannot penetrate with a 
power-law dependence to large distances but rather must 
decay exponentially with the decay length A c , because 
this component of the field induces Josephson currents 
to flow perpendicular to the layers. As the pancake vor- 
tex moves deeper into the stack, the Josephson coupling 
begins to play a greater role. When the pancake vortex 
is a distance Aj or greater from the surface, a Josephson 
core region of width 2Aj appears in the region between 
the vortex axis and the surface. Finally, at distances 
such that xo ^> A c , the current and field distribution can 
be characterized as a pancake vortex in which the fields 
at distances less than Aj from the axis are nearly the 
same as in the Josephson-decoupled case, and the mag- 
netic flux carried up through the pancake layer z = is 
4>q(s /2A\\). However, this magnetic flux does not flow ra- 
dially outward to infinity as in the Josephson-decoupled 
case but instead is confined within a highly elliptical 
field distribution consisting of an overlapping Josephson 
vortex-antivortex pair, which links the pancake vortex 
to the surface. Recall that when a straight vortex is at 
a distance x from the surface of an isotropic supercon- 
ductor of penetration depth A, the magnetic flux inside 
the superconductor, calculated accounting for the image 
vortex at — x, is <fio[l — exp(— x/X)]. As a pancake vor- 
tex moves from the surface to a position xq deep within 
the superconductor, it drags along a Josephson vortex 
(carrying magnetic flux in the +x direetion) whose axis 
is in the insulating layer at z — —s/2, and it also drags 
along a Josephson antivortex (carrying magnetic flux in 
the —x direetion) whose axis is in the insulating layer 
at z = +s/2. Accounting for the overlapping field dis- 
tributions, which nearly cancel each other, we find that 
the magnetic flux carried in the +x direetion through 
the space z < is </)q[1 — exp(— s/2A||)] « 0o(s/2A||); the 
same amount of magnetic flux is carried back in the —x 
direetion through the space z > 0. 

To give another example, consider an infinite stack of 
pancake vortices initially aligned along the z axis in an in- 
finite stack of Josephson-decoupled superconducting lay- 
ers. As diseussed at the end of Sec. II, the field and cur- 
rent distributions reduce to those of a line vortex in an 
isotropic superconductor of penetration depth A||4* The 
magnetic field is everywhere perpendicular to the lay- 
ers. Now imagine displacing all of the pancake vortices 
in the space z > s/2 by a distance xq in the x direetion, 
such that the pancake vortex stack now has a kink at 
z = s/2. In the absence of Josephson coupling, the re- 
sulting field and current distributions can be obtained by 
superposing those given in Sec. II. A component of the 



field parallel to the layers must arise in order to displace 
the magnetic flux </>o whose distribution is centered on the 
z axis for z <C — A| to a distribution centered on the line 
(x,y) = (a;o,0) for z ^> An. The component of the field 
parallel to the layers has a dipole-like distribution in any 
plane z = const, with a power-law dependence at large 
distances, but it decreases exponentially for \z\ > A| 
because of the sereening currents that flow parallel to 
the layers. In the presence of interlayer Josephson cou- 
pling, the above picture is altered, and it is now useful 
to think of kinked vortices as staeks of pancake vortices 
connected by Josephson strings (short pieces of Joseph- 
son vortices). The axes of the Josephson strings are 
confined to the insulating regions between superconduct- 
ing layers. As a consequence of the Josephson coupling, 
the radial component of the magnetic field is sereened 
on the length scale of A c by the induced Josephson cur- 
rents, which flow perpendicular to the layers. Although 
there is little perturbation of the original field distribu- 
tion when xq < Aj, the Josephson length, this situation 
is altered when xq > Aj, because in this case a nonlin- 
ear Josephson core appears along the string connecting 
the two pancake vortices centered at (x,y,z) = (0,0,0) 
and (x,y,z) = (x o ,0,s). The Josephson-energy cost of 
the Josephson string coupling the two semi-infinite staeks 
of pancake vortices is approximately (taking logarithmic 
terms to be of order unity) 1:L i 63 i 65 i 66 i 67 i 68 

E short (x ) « {c/)o/4ir) 2 xl/s\ 2 c ,XQ < \j = (X c /X ab )s, 

. (90) 

when the Josephson string is short and its core is not 
fully formed. The Josephson-energy cost is of the order 

fll. 61. 62. 63. 65. 66. 67. 68. 69 

Elong{x ) S» (<j) /4ir) 2 Xo/\ ab \ c ,X > Xj = (A c /A q6 )s, 

when the Josephson string is long and its core is more 
fully formed. However, it is not until xq ^ A c that 
a fully formed Josephson vortex (with width 2A C and 
height 2A a b) can streteh out between the upper and lower 
parts of the split stack of pancake vortices. In this case 
the energy cost of the Josephson string coupling the 
two semi-infinite staeks of pancake vortices reduces to 
(4>o H d, ab/ '47r) x , where 4>oH c i,ab/^ is the energy per 
unit length of an isolated Josephson vortex parallel to 
the layers and 61 i 67 i 69 i 70 

ff cl , a6 = — fV[ln(^)+1.55]. (92) 

is the lower critical field parallel to the layers. 

In anisotropic superconduetors consisting of 
Josephson-coupled superconducting layers, one may 
always regard the vortex structure as consisting of 
a superposition of 2D pancake vortices, which carry 
magnetic flux up through the layers, and Josephson 
vortices (or strings), which carry magnetic flux parallel 
to the layers but no net flux perpendicular to the layers. 
In transport experiments involving vortex motion, the 
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voltages are given by the Josephson relationsi^i The dc 
voltage parallel to the layers is Vu = (h/2e)v\\, where 
mi is the time-averaged rate with which 2D pancake 
vortices cross a line between the contacts, and the dc 
voltage perpendicular to the layers is proportional to 
V± = (h/2e)i>±, where u± is the time-averaged rate with 
which the axes of Josephson vortices (or strings) cross a 
line between the contacts. 

When the Josephson coupling is strong, vortex lines 
tilted with respect to the z (or c) axis can be described 
as tilted stacks of 2D pancake vortices or as a tilted lat- 
tice, where pancakes in adjacent layers are connected 
by Josephson strings. Such vortices, sometimes called 
kinked vortex lines^i*^ have been studied by numerous 
authorSi*2J2iiii2fai& However, when the Josephson cou- 
pling is very weak, a magnetic field applied at a small an- 
gle relative to the layers can produce a structure consist- 
ing of two perpendicular interpenetrating latticesS*2LZ& 
(called a combined lattice^ or crossing lattices^): a 
lattice of pancake vortices aligned nearly perpendicular 
to the layers and a lattice of Josephson vortices par- 
allel to the layers. The interaction between the two 
kinds of vortices leads to striking chain-like vortex pat- 
terns in highly anisotropic Bi-2212, which have been ob- 
served by Bitter decoration£2i2£ and scanning Hall-probe 
microscopyi^iSiiS Both techniques reveal the positions 
of 2D pancake vortices within about X a b of the surface. 
As shown by KoshelevjS in highly anisotropic layered 
superconductors the interactions between pancake vor- 
tices and Josephson vortices lead to deformations of both 
the pancake- vortex and Josephson- vortex crystals and to 
pinning of Josephson vortices by pancake vortices. 

At high temperatures and applied magnetic fields, the 
vortex lattice meltsriSiSiiSLSii and this process has even 
been directly visualizcd in Bi-2212 by scanning Hall- 
probe microscopyi 89 ' 90 The authors of Ref. |9(] used the 
formalism of Sec. IV A to infer the Lindemann parame- 
ter from the rms thermal fluctuations of pancake vortices 
vs magnetic field just below the melting transition. Much 
experimental and theoretical research has been devoted 
to vortex-lattice melting, and the reader is referred to 
reviews by Blatter et ali^I and Brandt 92 for a more com- 
plete discussion of this topic. 

The pinning of vortices by point defects is another 
topic where the interactions between pancake vortices 
and Josephson vortices play a key role. This difhcult 
subject is further complicated by the effects of thermal 
fluctuations, especially in the high-temperature super- 
conductors at the elevated temperatures where poten- 
tial applications are most interesting. The reader is re- 
ferred to the above reviewaSi^ and the recent paper by 
Kierfeld 9 ^ for further details about this subject. 



VII. SUMMARY 

In this paper I have presented solutions that permit 
the calculation of the magnetic-field and current-density 



distributions generated by a single 2D pancake vortex in 
an infinite stack (Sec. II), semi-infinite stack (Sec. III), or 
a finite-thickness stack (Sec. IV) of Josephson-decoupled 
superconducting layers. I have shown in Sec. V how to 
calculate the electromagnetic forces between two pancake 
vortices, and in Sec. VI I have diseussed some of the ways 
that interlayer Josephson coupling modifies the results. 

The results of this paper should be useful to those 
using probes (such as scanning Hall-probe microseopy, 
scanning SQUID microseopy, Bitter decoration, and 
magneto-opties) of the vortex-generated magnetic-field 
distributions above anisotropic high-temperature super- 
conductors. If the sample surface is parallel to the 
cuprate planes, these probes measure chiefly the mag- 
netic fields generated by pancake vortices within about 
A a b of the top surface. Although Josephson vortices (or 
strings) produce no net magnetic flux through the top 
surface, they can produce dipole-like stray fields if they 
are within X a b of the surface. On the other hand, if 
the sample surface is normal to the cuprate planes, such 
probes measure chiefly the magnetic fields generated by 
Josephson vortices within about A c of the sample surface, 
although pancake vortices within X a b of the surface can 
produce dipole-like stray fields outside the sampleiSi 

The pancake-vortex field and current distributions 
given in Sees. II-IV also could be useful in analyzing ex- 
periments such as Lorentz microscopy^S^SiSLSSiSSiiffi 
that probe the magnetic-field distribution throughout the 
sample thickness. 

Since the London model is at the heart of the above 
pancake-vortex caleulations, the resulting theoretical 
field and current distributions have unphysical singulari- 
ties at the pancake- vortex core, which is of size ~ £ a fc. 
Such singularities should have no experimental conse- 
quences for the above probes, which have insufficient res- 
olution to reveal details at this length scale. However, for 
probes of higher resolution it may be necessary to take 
into account the faet that the cireulating current density 
reaches a maximum at p w £ a {, and vanishes linearly as 
p — > 0, such that the singularity of the magnetic field 
at the pancake-vortex core is removed. The core effects 
could be treated approximately by using a vortex-core 
model that employs a variational core-radius parameter 
tv ~ Cab, as in Refs. I3H50H5I52I 
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APPENDIX: INTEGRALS USEFUL FOR D< A 

Scvcral intcgrals appear in Sec. IV B. Ali may bc cval- 
uated by starting f r0 mi2iii22 



du ~r, — = jPoW - y (z) 



(A.l) 



where H n (z) is the Struve function and Y n (z) is the 
Bessel function of the second kind (Weber's function), 
differentiating with respect to z, making use of recur- 
rence relations, integrating by parts, and making use of 
the properties tha1*i2i*i2S 



duJo(zu) 



duJi(zu) = — 



du 



J\ (zu) 



= 1. 



/o u 

The vector potential a^p, 0) is proportional to 



(A.2) 
(A.3) 



du 



J\ (zu) 



l + 1 -l[H 1 (z)-Y 1 (z)} (A.4) 



1,2 < 1, 

-,z> 1, 



(A.5) 
(A.6) 



where the limiting forms for z C 1 and z 3> 1 a re ob- 
tained from expansions givcn in Rcfs. 101 and 102. How- 
ever, Eq. 

may be obtained more simply by noting 
that, because of the properties of J\(uz), the integral 
when z <C 1 is dominated by values of u > 1, such that 



1 + u may be replaced by w; the resulting integral then 
takes the form of Eq. I|A.3|) . Similarly, Eq. (|A.6|I may be 
obtained by noting that when z> 1 the integral is domi- 
nated by values of u <C 1, such that l + u may be replaced 
by 1; the resulting integral may be evaluated using Eq. 
l|A.2jl . The limiting forms of the following integrals also 
may be obtained in a similar fashion. 

The magnetic field component b z (p, 0) is proportional 

to 



du 



uJo(zu) 



1 f 00 uJ\(zu) 

'Z l (T+~u)2 



(A.7) 



du (l-u)Mzu) 



1 



[H (z) - Y (z)} (A.9) 



1, 



(A.10) 



1 



,z > 1, 



(A.ll) 



and the net sheet current Ko(p) is proportional to 



du 



uJi(zu) 



2 
1 

z 

1 



du 0K ' A.12 
(1 + Uf 



[Htiz) -Y^z)] (A. 13) 



,z< 1, 



z > 1. 



(A.14) 
(A.15) 
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